Rejuvenation and Memory in model Spin Glasses in 3 and 4 dimensions. 
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We numerically study aging for the Edwards- Anderson Model in 3 and 4 dimensions using differ- 
ent temperature-change protocols. In _D = 3, time scales a thousand times larger than in previous 
work are reached with the SUE machine. Deviations from cumulative aging are observed in the non 
monotonic time behavior of the coherence length. Memory and rejuvenation effects are found in 
a temperature-cycle protocol, revealed by vanishing effective waiting times. Similar effects are re- 
ported for the D = 3 site-diluted ferromagnetic Ising model (without chaos). However, rejuvenation 
is reduced if off-equilibrium corrections to the fluctuation-dissipation theorem are considered. Mem- 
ory and rejuvenation are quantitatively describable in terms of the growth regime of the spin- glass 
coherence length. 

PACS numbers: 05.70.Ln,75.10.Nr,75.40.Mg 
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I. INTRODUCTION 

Nowadays the arena for comparisons between theory 
and experiments in spin-glasses physicsi is out-of equilib- 
rium dynamics^ Spin-glasses age^ as shown by the ther- 
moremanent magnetization: consider a spin-glass that 
has spent a time t^ below its glass temperature, T^, in 
the presence of a magnetic field. Let t be the time elapsed 
since the magnetic field was switched-off. The magneti- 
zation decays as a function of i/i^, even for i-u, ~ 1 day. 
The exponent ^ could be 1 {full- aging), "^ although there 
are recent experimental claims for ^ being smaller than 1 
{suhaging)fi' A somehow complementary experiment con- 
sists in keeping the system for tw below its glass temper- 
ature. Then, a magnetic field is switched on and the so- 
called zero-field cooled (ZFC) magnetization is recorded 
while it growsi^ Full aging can be also observed in the 
a.c. magnetic susceptibility, xi'^^ ^w) which, for a fixed to, 
decreases as tw grows. This time decay can be rescaled 
as a w-independent function of iot^,"^ although, experi- 
mentally, one is restricted to the range ujt^ > 1. 

Memory and rejuvenatior.yJ: are sophisticated manifes- 
tations of aging in experiments where the temperature is 
not kept constant. Rejuvenation arises when changing 
temperature from Ti to T2 (Ti and T2 smaller than the 
critical temperature, Tc) a system that has spent some 
time at Ti, so that xi^^tw) barely depends on tw Just 
after the Ti — > T2 change, aging restarts. The imagi- 
nary part of the susceptibility suddenly grows then re- 
lax. The iu, dependency oix"i^,tw) gets stronger, as for 
a younger system. Rejuvenation means that the relax- 
ation of x"('^7^u)) is very similar to the one of a system 
just quenched from T > Tc to T2. Sometimes it is said 
that the relaxation is identical to the one of a system 
instantaneously quenched to T2 from infinite tempera- 
ture (in Scct. lIIIBl below, we elaborate on the different 
meaning of instantaneous temperature quench in a experi- 
ment and in a computer simulation) . If the susceptibility 



just after the quench to T2 rises above the final value it 
had at Ti, one speaks of strong rejuvenation.^ On the 
other hand, when the system is put back at tempera- 
ture Ti, x"(w,tu,) continue its relaxation where it left 
it just before the temperature change {memory effect). 
These effects can also be observed in the real part of 
the susceptibility (see e.g. Fig. 1 of Ref. 0), although 
rejuvenation is very diminished as compared with the 
imaginary part. With the sophisticated dip- experiment 
temperature-change protocol,^ memory and rejuvenation 
are truly spectacular. 

Memory and rejuvenation have been found in sys- 
tems quite different from spin-glasses (see, however, 
Ref. llOll . Examples are structural glassesii, polymers 
(PMMAiii^) , and systems not particularly glassy (or 
not widely recognized as such), like colossal magne- 
toresistance oxidesii^ Moreover, a disordered ferromag- 
netic alloyfi^ becoming spin-glass at lower temperatures, 
has shown rejuvenation and memory, through the dip- 
experiment protocol (although in this case memory could 
be easily erased by lowering the temperature). Spin- 
glasses display the quantitatively stronger effects, but it 
is unlikely that the physical mechanism underlying mem- 
ory and rejuvenation are specific of spin-glasses. 

The above definitions for memory and rejuvena- 
tion need qualification. Under very small temperature 
< 5 X 10~^) the behavior of the 
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spin-glass is rather smooth. On the other hand, sharp 
memory and rejuvenation can be observed^'' for -^ — ~ 
0.07 . The crossover from small to drastic effects is ra- 
tionalized using effective isothermal waiting timesi^^ 
Consider the simplest temperature change protocol: a 
system is aged for time tw at temperature Ti, then its 
temperature is suddenly shifted from Ti to T2. After the 
shift, the ZFC magnetization is measured. The effective 



time, t! 



cff, shift 



is the age of the isothermally aged system 
at temperature T2, whose ZFC magnetization^ is most 
similar to the one of the temperature-shifted system (the 



two relaxations are not identicali^). Rejuvenation arises 
when i!j, '^ ' /t^j is below experimental resolution. 

Similarly, one can definai^ an effective time for the 
temperature cycle protocol Ti ^ T2 ^ Ti^ one keeps 
the system a time t^, at Ti , then shifts the temperature to 
T2, waits a time ^2 ~ 20t^, shifts back the temperature to 
Ti , switches on a magnetic field and then records the ZFC 
magnetization. The effective time t!p ''^^'^ ° is obtained by 
looking to the system aged at temperature Ti for a time 
tw + t'r '^^'^ ° whose ZFC magnetization is most similar 
to the one of the temperature-cycled system. One has 
memory, as we defined it above, when t^ ''^■^'^ °/t2 gets 
below experimental resolution. A large variety of spin- 
glasses experiments find-"^® for Ti > T2 



t: 



cff .cycle 



exp 



Ti~T2 



X0T2 



with^a 



a;o - 10" 



(1) 



(2) 



The theoretical investigation of these phenomena is 
less advanced than its experimental counterpart. Mem- 
ory and rejuvenation can be recovered in the dynam- 
ics of abstract energy-landscape models.^^ However, one 
wants to reproduce these phenomena in the Langevin dy- 
namics for the standard spin-glass model, the Edwards- 
Anderson (EA) model.— This dynamics for the EA 
model can only be investigated by Monte Carlo simu- 
lation. Yet, difficulties have arisen in numerical investi- 
gation of memory and rejuvenation&SL^SiSiSi Further- 
more, the progress achieved regards temperature-shift 
and temperature-cycle experiments. The dip-experiment 
protocol remains still as too complicated to be analyzed 
theoretically. 

Experiments where (T2 — Ti)/Ti is very small can be 
accounted for by the cumulative aging scenario piSiiS. con- 
sisting in the three following hypothesis: 

1. Aging is ruled by the growth of a coherence 
length, ^^ signaling the building of a spin-glass or- 
der. For isothermal aging, this length is named 
^T{t), t being the total time spent in the glass 
phase. This isothermal growth-law has been stud- 
ied in experiments^^ and simulations r^iSS although 
the measured ^rit) grows by an small factor in both 
cases. Numerically, a power law 
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fairly fits the data. However, more complicated 
rules have been usedi^*iSiiLi& 

2. The coherence-length always grows with time. It 
behaves continuously upon temperature changes. 

3. Effective times follow from the isothermal growth 
of the coherence length. Consider a temperature 



shift after aging for time t^ at Ti. One has 

CT,(i») = eT.(4f'"'). (4) 

A time t after the shift, the coherence length is 



r'"w = eT.(i+i' 
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Similar reasoning is used in the analysis of more 
complicated temperature-change protocols. 

Eq.Q is used in an indirect way, both in the analysis of 
simulations^:^ and experiments i^^i^^i^^ Relations such as 
(O, obtained in a different experiment, are used to con- 
vert the measured effective times into length-scales and 
viceversa. It is difficult to find in the literature direct 
data on the behavior of the coherence length upon tem- 
perature changes. A nice exception are the simulations of 
Ref . I21I where Eq.jSJ was directly checked. Those simula- 
tion spanned 10^ Monte Carlo steps (MCS). For compar- 
ison with experiments, recall that 1 MCS ^ 1 picosecond. 

Memory and rejuvenation appear as hardly compat- 
ible with the cumulative aging. Experiments showi^ 
that £,Ti{tw) < Ct2(^t2''' ' ) 'when the measured effective 
times are converted in length-scales, both for T2 > Ti 
and Ti < T2, in contradiction with Eq.Q). 

Two theoretical scenarios are currently being consid- 
ered to account for memory and rejuvenation. Rejuve- 
nation was interpreted in terms of temperature chaos^^ 
namely extreme sensitivity of equilibrium, states in the 
glass phase to small temperature changes. An overlap- 
length, lo{Ti,T2), is postulated to exist. Features at 
scale smaller than Iq are unaffected by a temperature 
change Ti -^ T2 while at larger scales the system is com- 
pletely reorganized. Rejuvenation is then attributed to 
large length scales and strong rejuvenation requires small 
Iq. The ghost-domain scenario (see Ref. llTl for a recent 
account) allows to reproduce memory in the chaos sce- 
nario. The other scenario^fi is closer in spirit to cumu- 
lative aging. Rejuvenation after a negative temperature 
shift would arise from the so-called fast modes involving 
length-scales smaller than ^Ti (tw), that were equilibrated 
at Ti but fall out of equilibrium at T2. Memory would 
arise from time and length scales separation: back to 
temperature Ti, fast modes re-equilibrate very fast so 
that aging continues from the previous Ti state. 

However, when it comes to actual calculations, it turns 
out that no convincing memory and rejuvenation has 
been found in computer simulations of 3D spin-glass 
models, either with a two-temperatures^iSLSSiSi or with 
a dip-experiment protocol. ^^ When the behavior of the 
coherence length is followed for times up to 10^ MCS,— 
Eqs.EI) and © are fulfilled even for ^^ ~ ±0.33. 
Consistently with this finding, when the temperature cy- 
cle protocol is analyzed in the EA modelrii2& the xo 
in Eq.(^) turns out to be of order 1 rather than of or- 
der 10~^ . Should xq not decrease significantly for larger 
times, the whole low-temperature phase of the EA model 
could be accounted for by cumulative aging (i.e. the low- 
temperature phase would not be a spin-glass phase). 



This contradiction with experiments is puzzhng. It 
could be indicating that the EA model lacks some crucial 
ingredient (maybe long-ranged dipolar interactions?). 
Or maybe memory and rejuvenation involve time and 
length scales unaccessible to present-day simulations. In- 
deed, experiments are performed on a time-scale which 
is about 10^ times longer than typical simulations. Yet, 
experimentally,^* there are around ^ 10^ spins in a co- 
herent cluster (hence £,T{tw) ^ 40 lattice spacings), while 
simulations achieve (see below) £,T{tw) ~ 10 lattice spac- 
ings. When length scales are confronted, the differences 
with experimental conditions do not seem so dramatic. 

As for higher space dimensions, a simulation^'^ of the 
temperature cycle protocol for the 4D EA model, yielded 
strong rejuvenation (as defined in Ref. ^). Yet, results 
in full agreement with cumulative aging, Eq.l^J, were 
reported for '^^^'^^ « ±0.125 (the simulation time was 
smaller than 10^ MCS). In the Migdal-Kadanof lattice ,i£ 
where rather larger times can be simulated, rejuvenation 
was found for ^i^^ ~ 0.1 , suggesting that a;o in Eq.l^l, 
does depends on the age of the system. 



In this work, we report simulations of a 3D (made with 
the SUE machine^) and a 4D EA model with binary 
(rather than gaussiar^^'^^'^^^) couplings. Our 3D sim- 
ulations are three orders of magnitude longer than previ- 
ous ones. We use real replicas^^ to study the coherence- 
length, that is directly calculated (not inferred from ef- 
fective times), through the temperature changes. In a 
temperature-cycle protocol, clear memory and rejuvena- 
tion effects are found for large values of '^''""^^ both in 
3D and in 4D f sections IIII Al and IIII B|l . We also observe 
strong rejuvenation (in the sense of Rcf. 8), but only 
if we neglect corrections to the Fluctuation-Dissipation 
theorem,^* The coherence-length is shown to decrease 
upon some temperature changes, in contradiction with 
cumulative aging. Moreover, a value t^ ''^^'^ '^ compatible 
with zero can be obtained for Ti — 0.9Tc and T2 — 0.4Tc, 
which is to be expected in view of Eqs.QJ and (|2Jl (sec- 
tions ^^21 a-nd ^^3 . Furthermore, we will show that 
the two-times dependency of the time correlation func- 
tion can be accounted for with surprising accuracy by 
the coherence-length at the two relevant times, both 
for isothermal aging and for temperature-shift protocols 
We perform exactly the same calcula- 



(section IIIID|) . 
tions for the 3D ferromagnetic site-diluted Ising model^ 
(where, in the absence of frustration, chaos is absent), ob- 
taining very similar results. Although temperature chaos 
is probably present in models for spin-glasses,^^'^^ our re- 
sults in the site-diluted ferromagnetic Ising model suggest 
that it plays no role in producing memory and rejuve- 
nation. This was maybe to be expected, since memory 
and rejuvenation is being found experimentally in mate- 
rials where chaos in temperature seems to be absent or 
where a thermodynamic glass transition has never been 
foundiiiii^ii^ Uufortunatcly, we have made no progressSS, 
in the analysis of the dip-experiment protocol. 



II. MODELS AND SIMULATIONS 

Specifically, wc consider Ising variables, ai — ±1, oc- 
cupying the nodes of a (hyper) cubic lattice in 3D and 4D 
with nearest-neighbor, quenched disordered interactions. 
We report results for the spin-glass with random ±1 cou- 
plings and the 3D site diluted Ising ferromagnet^- (spins 
are lacking with probability 1—p). We evolve the system 
using a sequential, local heat-bath dynamics. Our time- 
unit (1 MCS ~ 1 picosecond) is a full-lattice update. For 
the spin glass we studied the lattice size L = 60 in 3D 
(mostly in SUE), and L ^ 20 for 4D (on PC clusters) 
with some tests in L = 30 finding no differences. For the 
site-diluted model we studied L = 100. The number of 
disorder realizations vary within 16 and 240. 

In the following, we will call a direct-quench to the pro- 
cedure of placing a fully disordered system (infinite tem- 
perature) instantaneously at the working temperature. 
This corresponds to an infinite quenching-rate. 
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FIG. 1: Top: Naive x(^ = 27r/io,U), to = 5.3 x lO'' for the 
3D EA model vs. time. The T-cycle is Ti ^T2 ^Ti ^T2, 
each step lasting ts — 2 x 10*. The full line is a reference run 
at O.QTc- The inset shows the rejuvenation time (see text) vs. 
to. The dashed line is a direct-quench to T2. Bottom: as 



top part for D=4, to = 1.6 x 10 and ts 
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The Fluctuation-Dissipation Theorem (FDT) relates 
the autocorrelation function in zero magnetic field 

C{t^,tw + to) = y'^{<yi{t^)(Ti{t^ + to)) (6) 



to the real part of the susceptibility: xi^ — 27r/to, ^to) ~ 
[1 - C{t^„t^ + to,)]/T. Yet, off-equilibrium, FDT 
needs to be generalized replacing T by T/X[C] {X[C] 
is a smooth function^i of C{tw,tu, + to))- Hence, one 
assumes^iSSiSiSl to be in pseudo-equilibrium regime 
{bjtyj ^ 1 thus X[C] = 1), which is not always true. 
We also obtain spatial information from the correlation 



function of the overlap field, qi (t) = a. 



(1) 



W^f^W, 



built 



from two independently evolving systems with the same 



couplings, at the same temperature: 

C4(r,tu,) = —^{qi{tni)qt+r{tw)) ■ (7) 



III. RESULTS 
A. Strong rejuvenation? 

In Fig. ^ is shown the time-evolution of the naive 
X(w,t«,) (i.e. [1 - C{t^,t^ +tQ)]/T) for the EA model 
in 3D (top) and 4D (bottom) for a (double) temperature 
cycle: Ti ^ T2 ^ Ti ^ T2 (Ti = 0.9T, , T2 = 0.4rc). 
In 3D the system spends tg ^ 2 x 10'* MCS at each tem- 
perature (1000 times longer than previous works), while 
in 4D is = 10^ MCS. The results of a reference run, 
with temperature fixed to Ti, are also shown (continu- 
ous line). When the temperature drops to T2, x{LL>,tw) 
increases over the reference curve and starts a new relax- 
ation {strong rejuvenation^). When temperature is back 
to Ti, x(w,tu,) catches the reference run almost instan- 
taneously {memory). We call trej to the time that the 
rejuvenated x(w,i^) is above the reference run (see in- 
sets in Fig.^, that is found to grow consistently with to 
(much faster in 4D). It is then conceivable that an effect 
of macroscopic time-duration could be observed (experi- 
ments explore to ~ lO^'^ MCS). However, specially in 31?, 
^rcj < to. This implies that this strong rejuvenation is 
confined to the regime Lotw < 1, which is out of reach for 
measurements of the a.c. susceptibility (note that strong 
rejuvenation is not always observed experimentally in the 
real part of the susceptibilit'^) . 

In agreement with Ref. |23, the relaxing curve after the 
temperature drop is independent of ts on the explored 
range {ts = 10^ 2 x 10^ and 2 x 10** MCS in Z? = 3). 
Also shown in Fig. Q] is the relaxation of xi'^it^) for 
a direct-quench to T2 (dashed-line). Such an infinitely- 
fast temperature drop is not realistic (see section UlI B|l . 
Anyhow, the relaxation is not identical to the one after 
the temperature shift, but the two become very similar 
(in D — 3, this happens for i^, ^ 4io)- This is in marked 
contrast with previous simulations where tg ~ 10^ and 
to = 64. s For such a short times, one needs t^ ~ 500io 
for the two relaxation curves to approach each other. 

Yet, in Fig. ^ we assumed to be in pseudo-equilibrium 
regime. In order to estimate the FDT correction fac- 
tor X[C] we use the following procedure. We stay for 
time ts at Ti, then change temperature to T2, wait for 
time tiu and switch-on a small uniform magnetic field 
{h = 0.03). We then record the magnetization, m{tw+t), 
and C{tyj,tw + t). The sought X[C] factor is obtained 
drawing xT — m{t + tw)T/h versus C{tw,tw + t) (in- 
sets of Fig. 12)). The resulting plot, i^-independent for 
large t^r^ can be fitted with two straight- lines, yielding 
X[C]. In Fig. 121 (top) we show the time evolution of the 
x{ijJ,tw) for the same cycle as Fig. ^ with ts — 2 x 10^. 
Correcting with X[C] reduces rejuvenation to the point 
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FIG. 2: Top: As Fig.[T]for to = 4.3 x lO''^ and f « = 2 x 10^ but 
correcting FDT violations. Inset: Off-equilibrium fluctuation 
dissipation relation (see text). Data can be nicely fitted to 
two straight lines. Bottom: As in top part for T2 — 0.7Tc. 
Only the pseudo-equilibrium regime is explored in this case 
(inset). 



that strong rejuvenation is no longer seen. The suscep- 
tibility no more grows at the temperature drop to T2, 
although the relaxation restarts and still collapses ap- 
preciably with the x(a;, tw) curve obtained from a direct- 
quench. Similar conclusions are drawn in 4D.'^^ We re- 
port in Fig. 121 (bottom) results for a cycle with a smaller 
temperature step (Ti = 0.9Tc, T2 = O.TTc)- Here, (see 
inset), we stay in pseudo-equilibrium regime and reju- 
venation is stronger for the smaller temperature drop, 
once the correcting X[C] factor is considered. However 
the collapse with the direct-quench curve starts only for 
tw '^ 20to- 



1. The diluted ferromagnet 

Memory and rejuvenation have been found in other 
systems than spin-glasses. A disordered ferromagnetic 
alloy)^ becoming spin-glass at lower temperatures, has 
shown rejuvenation but much weaker memory, through 
the dip-experiment protocol. For comparison, we have 
simulated a site-diluted Ising model for p = 0.395 (Tc 
is accurately known^). All the interactions being fer- 
romagnetic, there is no temperature chaos in this sys- 
tem. We have simulated a. L = 100 system checking 
that S^T ^ L in our simulation window. We measure 
the naive susceptibility with the autocorrelation func- 
tion (O. Just for fun, we try a three steps protocol, 
Ti = 0.9T, ^ T2 = 0.7T, ^ Ta = 0.4rc ^ T2 ^ T3, 
staying ts — 10^ MCS at each temperature. The results 
are shown in Fig. (bottom) together with the results 
for an equal protocol for the 3D EA model (Fig. [31 top). 
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FIG. 3: Naive susceptibility in a 3 temperatures protocol (see 
text for details), for the binary EA D — 3 model (top) and 
the D = 3 site diluted ferromagnetic Ising model (bottom). 
The dashed line is a direct-quench to Ti. 
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In both cases rejuvenation and a double memory are ob- 
served. Also in two temperature cycles^ the susceptibil- 
ity behaves as in the EA model. To this level of anal- 
ysis, there is no clear difference bet"ween the Ed'wards- 
Anderson model and the site-diluted Ising model. 



B. Comparison with experimental direct-quench 

In view of Eqs.lQ) and Q, and the large temperature 
drop that we are studying, one would expect a perfect 
rejuvenation effect. However, Figs. ^ and [21 show that 
the relaxation after the first step at Q.9Tc considerably 
differs from the direct-quench (although this difference 
is smaller than for shorter simulationgSiSi) . This seems 
in plain contradiction with experiments (see e.g. Fig. 4 
of Ref . [l7|) . Yet, upon reflection, one realizes that the 
experimental direct-quench bears little resemblance with 
the simulational one. In fact, the experimental sample 
that is "instantaneously" quenched to 0.4 Tc, expends at 
least 10 seconds (~ 10^^ MCS!) in the spin-glass phase. 

In order to make a fair comparison with experiments, 
one should study the relaxation after a "soft" quench 
(Fig. 2)) from high-temperature to the working tempera- 
ture below the glass transition. Yet, the fastest quench- 
ing rate that can be achieved in experiments is far too 
slow to be reproduced in present-day computers. To 
achieve a very slow temperature drop from high tem- 
perature to working temperature, it is useful to consider 
Fig. n in a different way. One realizes that the system 
that has spent tg =2x 10* MCS at 0.9Tc, then suffers an 
instantaneous temperature drop to 0.4Tc is a better ap- 
proximation to the experimental direct-quench to 0.4rc . 
In fact, the system spends quite a long time close to the 
critical temperature, where the time evolution — recall 
Eq.(|3Jl — is faster. When looking to the double tempera- 
ture cycle in Fig.^ one needs to compare the relaxation 
in the first and in the second steps at 0.4Tc, the first 



FIG. 4: (Color online) Naive x(t^ = 27r/to,i™) of the 3D EA 
model vs. time, at temperature T2 — OATc, for several values 
of to and thermal histories. In all cases, tw is measured from 
the time of the (last) instantaneous quench to T2. Squares 
correspond to the instantaneous drop from infinite tempera- 
ture to T2. Circles correspond to the system that has spent 
ts = 2 X 10* MCS at Ti = O.gTc, then suffers an instanta- 
neous quench to T2. Crosses correspond to the system that 
has been a time ts at Ti, then time ta at T2, then time ts at 
Ti and finally suffers the instantaneous quench to T2. Dia- 
monds correspond to a gradual drop from 9Tc to T2 in 20000 
MCS (we incremented l/T in 0.113/7; every 10^ steps, the 
system spending 1.2 x 10* MCS in the spin-glass phase). 



corresponding to the reference direct-quench, the second 
being looked at as the temperature-cycled system. This 
comparison is shown in Fig. ^ together with the relax- 
ation after a soft-quench. 

The frequencies shown in Fig. 0] span three orders of 
magnitude. In all cases, the relaxation for the softly- 
quenched system^ that has spent 1.2 x 10'' MCS in the 
spin-glass phase, is much closer to the one of the cycled- 
system than the one of infinite quenching rate. Further- 
more, the relaxations for the first and the second steps 
at 0.4Tc are identical, up to our statistical accuracy (see 
Fig. [SJ. This you may wish to call perfect rejuvenation. 

In Fig. [51 we perform a detailed comparison between 
the soft-quench (with two quenching rates) and the two- 
steps protocol. To have a feeling of the frequency de- 
pendence, we show the smallest and the largest frequen- 
cies in Fig. 21 For very short times, in the two-steps 
protocols we find a quick decay of the susceptibility due 
to the sharp temperature drop. On the other hand, the 
softly-quenched system shows a basically constant behav- 
ior (the slower the quench, the lower the intial plateau 
is). When time becomes of the order of the total time 
spent in the spin-glass phase during the soft-quench, the 
susceptibility starts to decay and becomes very similar 
to the two-steps protocol. At to = 8192, the two soft- 
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FIG. 5: (Color online) Close-up to the top and bottom panels 
of Fig. 2] excluding the data for the instantaneous drop from 
infinite temperature to T2 = OAT^. We also include data 
(triangles) for a still slower drop from 9Tc to T2 in 200000 
MCS (we incremented l/T in O.llS/Tc every IC* steps, so that 
the system spends 1.2 x lO'' MCS in the spin-glass phase). 



quenches catch the relaxation of the two-step protocol 
and become identical. At the smallest frequency, the 
fastest quench approaches but does not catch the two- 
steps relaxation. On the other hand, for the smallest 
quenching rate, the relaxation curve becomes identical 
to the one of the two-steps protocol for tw ^ ^Oi which 
corresponds to the experimentally accessible time range. 
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FIG. 6: (Color online) As in Fig.|ll for the diluted ferromag- 
net and tg = 10^ MCS. Diamonds correspond to a gradual 
drop from 9Tc to T2 in 20000 MCS (we incremented l/T in 
O.m/Tc every 10^ steps). 



Quite similar results are obtained for the diluted fer- 
romagnet, as we show in Fig. |^ Although it cannot be 
noticed at the scale of this plot, the relaxations in the 
first and second temperature-step are not identical for 
the Ising model. They may be made to collapse if the 
times in the second step are rescaled by a factor of 1.2 (in 
a protocol O.TTc -^ OAT^ -> 0.7Tc -^ OATc, the needed 
rescaling factor is 2). 



C. The coherence-length 

The coherence-length may play a crucial role^° in this 
physics, and should be followed in detail during temper- 
ature changes. This was done previously in Ref. l2ll for 
times up to 10^ MCS. Results in agreement with Eq.® 
were reported. We show here qualitatively different re- 
sults for our longer simulations in 3D. 

The coherence-length may be obtained from non 
self- averaging integrals of C4{r,tw) using a second- 
momentum estimator's*^ Not having so many samples 
at our disposal, we have obtained C4(r, t^) (which is self- 
averaging for not very large r). The resulting curve has 
been fitted toS^ 



C4(r, tu,) = — exp 



atn 



(8) 



In 3D, we find fair fits in the range 2 < r < 20, fixing 
a = 0.65 and (3 = 1.7 for all times and temperatures. The 
constant behavior of a does not agree with the results for 
the 4D model with Gaussian couplings,S^ To estimate 
errors in the three parameters fit JSJ is very difficult. 
To have a feeling of their magnitude, let us report that 
a = 0.7 yields good fits as well, with a 10% increased £, 
estimate. 

See in Fig. |3 (top), S,{tw) for a direct-quench to T2 = 
OATc and for a thermal cycle Ti ^ T2 ^ Ti with Ti = 
0.9Tc and tg — 2 x 10^. A power law with exponent 
^ 0.144 fits nicely £,Ti (tw) for t,i, < t^, while the exponent 
for the direct-quench to T2 is r^ 0.065 (full-lines in Fig.[7J- 
top). Note that the exponent follow Eq.©. During the 
T2-step, ^ grows over the Ti value, and it is larger than 
for the direct-quench to T2. However, ^ decreases when 
the system is back to Ti . Memory is striking: data for 
the second Ti step, if translated back ts MCS, are on top 
of the fit (obtained for t^ < tg !). Let us stress two points 
regarding this result: 

1. The coherence-length can decrease upon temper- 
ature changes, violating cumulative-aging, Eq.©, 
and in contradiction with the time and length 
scales separation scenarioi22iS However, the effect 
is not symmetrical for negative and positive tem- 
perature shifts, as it was inferred experimentally 
from effective-time measurements^ 

2. The effective time for the temperature cycle is com- 
patible with zero (within our accuracy). This im- 
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FIG. 7: Coherence-length vs. time, for the 3D EA 
model (top) and the diluted- ferromagnet (bottom). The 
thermal history has been a T-cycle Ti —> T2 ^ Ti, being 
ts = 2 X lO'' for the EA model and ts = 10^ for the Ising- 
diluted model. (Crosses; second Ti-step data, translated back 
in time ts ; Dotted lines: fits to ^(i^) — At^ for i„ <ts). 



FIG. 8: Top: C(i™,i„ +to)(C(io + t™)/5(t™))^/^ (points) and 
C{tw,tm +to) (lines) vs. tw, for the direct-quench to T2, for 
the EA model in 3D. Inset: As main plot for the T2 step of 
the T-cycle of FigQ Bottom: As inset of top panel for the 
diluted Ising model. 



plies that, for tg ~ 10^, x^ in Eq.(^ is not of order 
one, as it was foundSi for ts ~ 10'' . 

As before (bottom part of Fig. [Tjl, the behavior of the 
diluted ferromagnet is completely analogous to the one of 
the EA model (including the power-law growth of ^t(O)- 



D. The coherence length and two times 
correlations 



A rather crucial feature of aging^ is that two time- 
scales, io and tw, are involved. One would like to relate 
the one time quantity £,Titw), to the two times correlation 
function. A crude estimate for io 3> t^ is 



Cit^^ tw -f io) oc 



e^/2(i.+to) 



(9) 



i.e. the coherent cluster that at time t^ + io has linear 
size £^{tw + to), at time t^ was composed of mutually 
incoherent clusters of linear size C(^t«) • 

Indeed, (Fig.il top), the factor ^^/^(t^ + to)/^^^'^{t^) 
absorbs almost all the t^ and to dependency of C{ttu, tw + 
to), both for a direct-quench to T2 and for the T2 part of 
the thermal cycle. Note that even the constant value for 
C(iu,,t^-t-to)C^^^(^w)+io)/C^^^(i«;) is equal for the direct- 
quench and for the thermal cycle. Also at Ti, Cit^, t^ + 
to)^^^^{tw + to)/^^^'^{tn,) is constant within a band of 
width 5% of its mean-valuei^ Quite similar results are 
obtained for the diluted ferromagnet, as we show in the 
bottom part of Fig. |H1 In spite of the crudeness of the 
argument leading to Eq.Q and the uncertainty in the 
determination of ^, the results are surprisingly clear. 



Note that if Eq.lO was exact the dynamics would be 
of the one-sector typep^ which we do not believe to be the 
casei^ Anyhow, given Eq.®, full-aging^ is natural for a 
power-law growth of ^(i^,). 

Thus, memory and rejuvenation are driven by the rate- 
growth of ^T(iio) rather than by its value or by the short- 
distance behavior of C4(r, tu,).^°'^^ In our simulation, re- 
juvenation is due to a growth of ^ upon cooling (probably, 
because of a sudden fall into a nearby energy minima), 
provoking a change in the evolution of C(tw,tw + to). 
When temperature is shifted back to Ti, £,ti continues its 
growth as if it had never being at T2 with analogous con- 
sequences for the correlation- function (memory). This 
implies a non monotonic behavior of S.T{t), in contradic- 
tion with cumulative-aging, Eqs.Q and jSJ. 



E. Results for tutn, < 1 

Measurements of a.c. susceptibility are usually con- 
fined to the region ujtm > 1. On the other hand, thermal 
magnetoresistance measurements can yield information 
on the time regime to ^ t^^ . It is therefore worthwhile to 
have a look to our correlation functions in this regime. 

The main issue here is the characterization of the time 
decay of correlations (see Ref. |33 for a recent study). 
One finda^ that, at least when the correlation function 
lies in some intervals (say Ci < C{tw,tw + to) < C2), 
it behaves as a function of to/tf^. It is possible that 
diiferent intervals for the correlation functions (usually 
called time-sector^) are ruled by different exponents. 
The existence of more than one time-sector is a nec- 
essary (but not sufficient) requirement for dynamic 



ultrametricityJ. Some indications of the presence of more 
than one time-sector in the dynamics of the EA model 
in 3D were found in Ref . ISft 
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FIG. 9: Correlation function C(tu,,iu, + io) scaled by ijj"'^ (see 
Eg. 1101 ') after a direct-quench to T = O.STc, versus to/t^. 
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FIG. 10: (Color online) Correlation function C{tw,tw + io) 
after a direct-quench to T = GATc, versus io/iiu- Inset: Data 
of main plot, versus t/tJJ,, /i = 2/3. We also plot data for a 
system that has been 2 x 10* MCS at O.OTc, then suffers a 
temperature shift to 0.4rc . 



is controversial. 
= 1 (full-aging), 
were ascribed to 
In the 



The experimental value of fx 
Recently^ it was claimed that it is /z 
Previous failures in recognizing this, 
quenching-rate effects^ (recall also Sect. IIIIBfl . 
quickest possible quench, ^ = 1 was clearly identified. 
This interpretation was recently disputed by the Saclay 
group, ^ that find /^ < 1 . 

Regarding computer simulations, the following scaling 
form for the time correlation function was proposed;^ 
and found to work for a restricted ^o f^nd t^ range: 



C(t^,,<^-hio) = to'''^'* 



(10) 



This equation implies that full aging should be observed 
in the io ^ iw regime. More recent and longer simula- 
tions'^^ found deviations from Ea. (|l()|l . at least at some 
temperatures. 

In Fig. 121 we plot our data for a direct-quench from 
infinite temperature to O.STc. Ea. ljlOII works nicely with 
x{T) — 0.02, which can be interpreted as evidence for 
full-aging behavior. However, at T = QATc (see Fig. llO|l . 
we have being unable of finding a working x{T) . Ac- 
tually, the relaxation is better interpreted in terms of 
two time-sectors: for to ^ tw, ^ = ^ seems to provide a 
proper scaling, while, at shorter io, M = 2/3 does a better 
job. 

In the inset of Fig. ^| we compare the relaxations 
aX T — QATc for the direct-quench and for the sys- 
tem that stayed 2 x 10* MCS at 0.9Tc, then suflFers a 
temperature shift to 0.4Tc (recall that this is our slow- 
est quenching-rate from infinite temperature). For the 
slowly quenched system, the relaxation belongs to the 
/I = 2/3 time-sector. Moreover, the relaxation after the 



T-shift provides a limiting curve, in the large tw limit, 
for the /I — 2/3 piece of the direct-quench relaxation. 

At least within the time range that we can study, it 
seems that the quenching-rate does exert influence in the 
measured value of /i (as proposed in Ref. 0). 



IV. CONCLUSIONS 

In this work, we have compared memory and rejuve- 
nation effects for the Z? = 3 and D = A binary EA model 
and for the D — 3 the site-diluted Ising model, finding 
quite similar results. Maybe the most important ques- 
tion we have addresed is whether the Edwards- Anderson 
model has a spin-glass low temperature phase, or not. 
The longer time scale reachable with the SUE machine^ 
(three orders of magnitude longer than previous work), 
has allowed us to conclude that spin-glass behave like 
experimental spin-glasses in a number of ways: 

1. The relaxation of the a.c. susceptibility after a 
large temperature shift is as for a direct-quench, 
provide that one does not interpret the term direct- 
quench literally in the simulations (Sect lIIIB|l . We 
find little dependency in the previous thermal 
history, once microscopically short times are ne- 
glected. 

2. Values of a;o, defined in Eq.(^), are not of order 1 for 
waitin g tim es ~ 10^ MCS (contrary to the finding 
of Ref/lH for waiting times of order -- 10^ MCS). 
This time dependence of xq is consistent with the 
findings in the Migdal-Kadanof lattice.^ 



3. The coherence length may decrease upon temper- 
ature changes, in disagreement with the cumula- 
tive aging scenario, and in agreement with recent 
experimentsiiSiii 

4. The growth-rate of the coherence length rules the 
decay of the time-correlation functioni^ A very 
simple formula, Eq. @ , accounts for the behavior of 
C(tyj^ ty^+to) with surprising accuracy, for different 
temperature protocols. Contrary to the findings of 
short 4D simulations)^^ we do not find a strong tem- 
perature dependency of the replica-field correlation 
function, Eq.Q, at short distances. We rather as- 
cribe rejuvenation to the coherence-length, which 
rules the long-distance behavior of the correlation- 
function. 

5. The exponent /i for the io/^(i scaling of the decay of 
C{tu],tw + to) , depends on the quenching-rate. In 
agreement with recent experiments^ we find that 
the faster the quench is, the easier it becomes to 
obtain /i = 1 (see, however, Ref. l3 for experiments 
contradicting this expectation). 

The above findings suggest that the 3D EA model be- 
haves as experimental spin-glasses do, contrary to what 
was inferred from previous shorter simulations ASLSSi^ 
However, there are a few important differences. First, 
the out-of-phase and the in-phase susceptibility behaves 
quite differently in experimental spin-glasses (see e.g. 
Ref 1^ . This seems not to be the case for the EA 
modelfSiS^iS within the accessible time-window. Sec- 
ond, the coherence-length is not found to decrease un- 
der positive temperature shifts, in contradiction with 
the inferred behavior from experimental measurements of 
effective-times ^Siii Third, in the dip-experiment protocol 
we have found^S no memory and very weak rejuvenation, 
in agreement with Ref. 22 even if we study frequencies 15 
times smaller. Fourth, recent simulations do not^^'^^ find 
stronger memory and rejuvenation effects for Heisenberg 
than for Ising models of spin-glasses. This finding is in 



contradiction both with experiments (see e.g. Ref lTq) and 
with a recent calculation of overlap lengths in the Migdal- 
Kadannoff approximation^^ The physical reasons under- 
lying these differences between our best model for spin- 
glasses and experimental spin-glasses are not yet under- 
stood. It is of course possible that longer times need to 
be studied. However, when time scales are converted to 
length-scales, one finds that numerical coherence-lengths 
are smaller than the experimental ones only by a small 
factor (see also Ref. l43|) . 

We have shown that the ferromagnetic site-diluted 
Ising model (where chaos is absent) follows very closely 
the behavior of the EA model, at least in the ujtw > 1 
regime. This is not totally unexpected, as we already 
know experimentally that systems quite different from 
spin-glasses show memory and rejuvenation. ^^'^^■^'^•^"^'^» 
The natural conclusion is that chaos, although probably 
present in realistic spin-glassesjSLSii^ need not be in- 
voked to explain memory and rejuvenation. However, it 
has been arguediS that the experimental protocol intro- 
duced in Ref. Il6| may help to discriminate temperature 
chaos from a somehow trivial restarting of the dynamics 
(yet, see Ref. 14^ . More work will be needed to assess the 
usefulness of this new classificationr^ of aging systems. 
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